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1. Introduction 
Let (X, Y) be a bivariate random vector with stochastic representation 
(1) {X,Y)^{RUi,RU2), 



Abstract. Let (X, Y) = {RUi, RU2) be a given bivariate scale mixture random vector, with R > being independent 
of the bivariate random vector {Ui,U2)- In this paper we derive exact asymptotic expansions of the joint survivor 
probabihty of {X, Y) assuming that R has distribution function in the Gumbel max-domain of attraction and (f/i, U2) 
has a specific tail behaviour around some absorbing point. We apply our results to investigate the asymptotic behaviour 
, of joint conditional excess distribution and the asymptotic independence for two models of bivariate scale mixture 

^Nj ' distributions. Furthermore for our models we derive an expression of the residual dependence index rj. 

u 

< 

Oh 

. where i? > is independent of the bivariate random vector (C/i, U2) (= stands for cquahty of distribution fmictions). 
^ ' The random vector {X, Y) has a scale mixture distribution; a canonical example of such {X, Y) is a bivariate spherical 
random vector with rotational invariant distribution function with (C/i, J72) uniformly distributed on the unit circle 
of K^. In this model (see Cambanis et al. (1981)) the dependence between Ui and U2 is a functional one, namely 
\^ ' Ul + U2 — 1 almost surely, and 
^ '. 

0\- (2) (C/i,[/2) = ihWj2\/l-W^), 

o 

Tjj- . with W e (0,l),/i,/2 e { — 1,1} almost surely, beta distributed with parameters 1/2,1/2 and P{Ii = 1} = 
P{l2 = 1} = 1/2. Furthermore, Ii, I2, R,W are mutually independent. 

If R is such that R^ is Chi-squared distributed with two degrees of freedom, then X and Y are independent Gaussian 
^ random variables with mean zero and variance 1. 

: Our main interest in this paper is the tail asymptotics of the joint survivor function of {X, Y). For Gaussian random 
^ ' vectors the asymptotics of the joint survivor probability is well-known, see e.g., Berman (1962), Dai and Mukherjea 
■ " " ' (2001), Hashorva (2005), or Lu and Li (2009). Results for elliptical and Dirichlet random vectors are obtained in 
Hashorva (2007, 2008, 2009c) and Manner and Segers (2009). Note that the elliptical model is derived by extending 
@ to 



(3) (C/i,t/2) = ihWJipW + hp.Vl-W^), pe(-l,l), p.:=V^ 



where the additional parameter p corresponds to correlation coefhcient of X and Y if R^ is Chi-squared distributed. 
Hashorva (2007) generalises the known asymptotic results for Gaussian random vectors to the more general class of 
elliptical ones by exploiting the fact that the asymptotics of the joint survivor probability is primarily determined by 
the asymptotic properties of the survivor function _F 1 — of the associated random radius R. 
Specifically, in the aforementioned paper the principal asymptotic assumption is that F is in the Gumbel max-domain 
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of attraction, which means that for some positive scahng function w 

F(x + t/w(x)) , , 

(4) Urn ^ —/ ^ ^ " = exp(-i), e M. 

x^oo F{x) 

As shown in Hashorva (2007) condition (j4]) is crucial when {X, Y) is an elhptical random vector with stochastic 
representation ([T]) . More specificahy, by the aforementioned paper for any a G (p, 1] 

(5) P{X > x,Y > ax} - — -—^—-F{apx), 



27r(l — ap){a — p) v{apx) 



where 



(6) flp := p^^ \J \ — 2ap -\- > 1, v{x) :— xw{x), x £ R. 

Throughout this paper f{x) ^ g{x) means Yanx^oo f{x)/g{x) = 1, and (|4]) is abbreviated by F G GMDA{w). 

If is beta distributed with positive parameters a, /3 (its distribution function is denoted by beta{a, /3)), then {X, Y) 
is a generahsed Dirichlet random vector. Hashorva (2009c) extends ^ for the class of Dirichlet random vectors. As 
indicated in Hashorva (2009b) for certain asymptotic problems the distributional properties of {Ui, U2) do not need 
to be explicitly known. A natural question that arises concerning the asymptotics of the joint survival probability of 
X, Y is that if F G GMDA{w) what models for the dependence between ([/i, U2) would lead us to asymptotic results 
similar to ([5])? 

In this paper we answer the above question for two specific models: The first one is refereed to as the unconstrained 
dependence model, or simply Model A. In that model we assume that Ui G (0, 1] almost surely, and further impose 
an asymptotic assumption on the behaviour of the {Ui,U2) around some absorbing point (l,a) (see ([9|) below). 

The second model (or simply Model B) motivated by ^ is referred to as the functional dependence model. More 
specifically we assume the stochastic representation 

(7) (f/i,C/2) = ihW,phW + l2z*iW)), pG(-l,l), 

with z* some positive measurable function, W G (0, l),/i,/2 G { — 1,1} almost surely, and Ii,l2,W are mutually 
independent. 

We present three applications of our results: The first one establishes an asymptotic approximation of the joint 
conditional excess distribution. In the second application we discuss the Gumbel max-domain of attraction of bivariate 
distributions related to our Model B. In our last application we derive an explicit expression of the residual tail 
dependence index rj for bivariate scale mixture random vectors extending a recent result of Hashorva (2010) for 
elliptical random vector. 

Organisation of the paper: In the next section we state our first result dealing with some general scale mixture 
bivariate random vectors which fall under Model A. We introduce in Section 3 some constrains on the dependence 
function of (C/i, U2) via ([7]), and then investigate the tail asymptotics of interest for Model B showing a generalisation 
of ^ in Proposition [21 Three applications of our results are given in Section 5. Proofs of all the results are relegated 
to Section 6. 



2. Tail Asymptotics Under Unconstrained Dependence 

Consider a bivariate scale mixture random vector {X, Y) — {RUi, RU2), where R has distribution function F (denote 
this i? ~ F) satisfying (|4]) with some positive scaling function w. We assume throughout this paper that F has an 
infinite upper endpoint. Hence by dU (see e.g., Resnick (2008)) 

(8) v{x) :— xw{x) 00, a; — ?► 00. 

Given a constant a G (0, 1] we investigate the asymptotics of 

Pa,s,-n-x P{X > x[l + S/v{x)],Y > ax[l + ri/v{x)]}, a; -)■ 00, 
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for any (5, 77 £ [0,oo). The re8-son for dcs-ling with the 8-syniptotics of Pa.s^ri^x 

is our interest concerning the approxi- 
mation and estimation of the joint conditional excess distribution, see the first apphcation in Section 4. 
Throughout in the sequel we assume that Ui is a bounded random variable. Without loss of generality we consider 
only the case Ui has distribution function with upper endpoint equal 1. This implies that Pa.s,ri;x < for ^ 

positive. For both Model A and B we show below that this upper bound is too crude; roughly speaking we have the 
asymptotic behaviour 

Pa,5,r,;x ^ 'lJj{x)F{x), 

with ip some positive function decaying polynomially fast to as a; — 00. 

In addition to the Gumbel max-domain of attraction assumption, we impose next a certain asymptotic behaviour of 
{111,1/2) around (l,a), namely 

/o^ r f {£^1 > 1 - (.s - S)/x, U2 > a(l - {s - r))/x)} . f x \ ^ rn ^ w ^ ^ 

(9 hm — -j = Cais,S,T]), V(5,ry e 0,00 ,Vs e 0,00), 

"{Ua > 1 — ^/x\ 

with a positive measurable function and Ua ■= min(C/i, 1/2/0.). If 5 > 77, then 

^a{s,6,T]) ^ £,a{s,6 ~r],0), £,a{s,T],S) = ^a{s,0,S -7]), Vse(0,Oo). 

Further, for (5 = r/ = 

(10) P{Ua>l-s} = s''La{s), Vs>0 

holds for some 7 S [0, 00), with La a positive measurable function such that \imsio La{s) / La{ts) — l,Vt > 0, i.e.. 
La is a slowly varying function; see Bingham et al. (1987), Embrechts et al. (1997), Falk et al. (2004), De Haan and 
Ferreira (2006), Jessen and Mikosch (2006), Resnick (2008), or Omey and Segers (2009) for more details on regularly 
varying functions. 

Next, we formulate our first result. 

Proposition 1. Let (X,Y) = (RUi, RU2) be a bivariate scale mixture random vector with F a positive random 
variable being independent of {Ui, 1/2)- Suppose that F has an infinite upper endpoint satisfying (|4]) with some positive 
scaling function w, and Ui G (0, 1] has distribution function with upper endpoint L If a ^ (0, 1] is such that ([9]) holds, 
then for any 5,rj £ [0, 00) we have 

"Fix) 

(11) Pa,S,ri-x ^ Js.,r,La{l/v{x))j^^j^-^, 

with La satisfying pU|) and 

/"OO 

A') / Ca(s, 5, 77) exp(-s) ds e (0, 00). 



Remarks: (a) In view of Lemma 6.1 in Hashorva (2009b) for any A S (l,oo),c G M and F as in Proposition [T] we 
have 

(12) hm = 0. 

x^oo F{x) 

In fact (IT2|l follows directly from Proposition 1.1 in Davis and Resnick (1988), see also Embrechts et al. (1997) p. 586, 
and Al. in Hashorva (2009c). 

Further we have the self- neglecting property of w, i.e., 

(13) w{x + t/w{x)) ^ ^ 

w(x) 

holds locally uniformly for t G M. Refer to Galambos (1987), Reiss (1989), Embrechts et al. (1997), Falk et al. (2004), 
de Haan and Ferreira (2006), or Resnick (2008) for details on the Gumbel max-domain of attraction, 
(b) Under the assumptions of Proposition [T] it follows that 

(14) Hix), 
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where H is the distribution function of RWa with Wa a positive random variable independent of R such that 

P{Wa>i-s}^JsMs)s''^ se(0,i). 

See also Example 1 below. Furthermore (fTTj) holds locally uniformly in 77, S. 

(c) Since for S — r] — ([T0| holds, then Jo,o = r(7 + !)■ By the monotonicity of Js^rj in S, rj we obtain 

J5,r, <r(7 + l), V5,?7>0. 

(d) By dH), ((ni) and ((I3|) it follows that 

Js,r, = exp(-J) Jo,,,_5, V77 e [0, (5], 
which follows also directly by the definition on Jg^t] and 
We present next three illustrating examples. 

Example 1. Let Ui,U2 be two random variables taking values in [0, 1] such that U2 > Ui almost surely. Suppose 
that P{Ui > I — s} ~ s''L{s), s e (0, 1) with 7 e [0, 00) and L a slowly varying function at 0. Since for any a; > 1 

P{Ui > 1 - l/x, U2> l-l/x} = P{Ui > 1 - l/x} 

if R is independent of {Ui, U2) satisfying the assumptions of Proposition [1] we obtain 

(15) P{RUi > x} = P{RUi> x,RU2> x} - r(7 + 

We note that for Ui ~ beta{a, (3) the asymptotics in ([TSl) is shown in Berman (1983), see also Berman (1992). For 
the more general case that Ui has a regularly varying survivor function see Theorem 3.1 in Hashorva et al. (2009). 

Example 2. (Linear Combinations) Let Si ~ Gi,i = 1,2 be two independent random variables with values in [0,1] 
such that 

(16) limg^T^iM = V.>0,z = l,2, 

2:^00 Gt{l - l/x) 

with 7i G [0, 00). Let Ai, A2 G (0, 1) be given constants with Ai > A2, and set 

Ui := XiSi + X.1S2, Xi'.^l-Xi, i = l,2. 

By the definition both Ui,U2 have upper endpoint 1. For any ^, 77 > we have (the proof is postponed to the last 
section) 

2 

(17) P{Ui>l-{s-5)/x,U2>l-{s-ri)/x} - f(s, ,5, 77) J] G.(l - l/x), s>max(5,?7), 

i=l 

with ^ 

|(s,5, 77) := 72 j (^max(^0,min([s - (5 - Aiz]/Ai, [s - 77 - A2z]/A2)j j z^^~^ dz. 

Note that for s G (0, max((5, 7;)] p7|) holds with ^(s, 5, rf) — 0, and when 5 = -q — Q 

(18) 1(5,0,0) ~ C^ui.MMs'''^^\ 
with 6*71,72, Ai.Aa e (0,00) given by 

C^^..l^MM -V / [1 - ATtp '^^ + A2"^' / ' \l-X'2tVH^^-^dt. 
Jo Jl 

Consequently, ^ holds with 

^i(s,(5,77) := 6*71,72, Ai,A2'f(s, 5, 77)s"''i"'^n(^>„iax(5,,,)), s > 0, 
where 1() is the indicator function. Thus with R~ F such that F G GMDv4(i(;) the result of Proposition [1] holds. 
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Example 3. (Farlie-Gumbel-Morgenstern Dependence) Let Ui ^ Gi,i — 1,2 he two random variables with values in 
[0, 1]. Suppose that for some KG [0, 1) 

P{Ui>x,U2>y}^Giix)G2{y)[l + KGi{x)G2iy)], Va;,y e [0,1]. 

The bivariate random vector {Ui, U2) possesses thus the Farlie-Gumbel-Morgenstern distribution (see for more details 
Hashorva and Hiisler (1999)). If p6|) holds, then for any (5, 77, s e (0, 00) we obtain 

2 

P{Ui > 1 - (s - 5)/x, C/2 > 1 - (5 - fi)/x} ^ (s - S)l'{s - r,)l' J]G,(1 - 1/x), 

i=l 

with (a;)+ := max(x, 0), a; G M. Consequently, if the positive random variable R F is independent of (Ui, U2) and 
F G GMDA{w), then locally uniformly in 6, rj 

POO 

(19) PaAv;- ~ / (t-(5)^(t-?7)Texp(-t)di[]G,(l-l/«(a;))F(x). 

•^0 i=i 

For any a G (0, 1) we observe another asymptotic behaviour, namely 

P{Ui > 1 - (s - 5)/x, U2 > a(l - (s - v)/^)} -{s- S)l'Giil - l/x)G2(a). 

Consequently, 

~ G2(a)r(7i -I- 1) exp(-(5)Gi(l - l/w(a;))F(a;). 

3. Tail Asymptotics For Functional Dependence 

In this section we deal with bivariate scale mixture random vectors assuming that the dependence between the 
components is determined by some deterministic function. Explicitly, let (X, y) be a bivariate random vector with 
stochastic representation 

(20) {X,Y) = {RhW,pRhW + Rl2Z*{Wj), pG(-l,l), 

with (/i, 12), -R > 0, W G (0, 1) mutually independent, and z* : [0, 1] [0, 1] a positive measurable function. 
We assume that the distribution function F of R has an infinite upper point, and that of W has upper endpoint equal 
1. In the sequel /i,/2 take values in { — 1, 1} with P{Ii = /2 = 1} £ (0, 1]. We allow Ii and I2 to be independent. 
The random vector {X, Y) is a scale mixture random vector for which the dependence of the components is being 
determined by p, z* and the random variables R, W, Ii, i — 1,2. We refer to the implied dependence of the components 
as the functional dependence. Note in passing that if 

VF^ ~ 6eto(l/2,l/2), z*{x) = ^/l-x^, a;G[0,l], 

and /i,/2 are independent assuming values —1, 1 with equal probability, then {X,Y) is an elliptical random vector. 

Generally speaking, under the setup of (|20|) it turns out that the local asymptotics of the density function of W is 
important. More precisely, we are able to provide an asymptotic expansion of 

Pa-x ■= Pa,o,0:x = P{X > x,Y > ax} , a G (0, 1], a; > 

requiring further that 

(21) P{W^l/ape{Kiu,K2u)} ^ Lk,,kAuW% 7ae[0,oo) 

holds for all u > small with Ki < K2,Ki,K2 G M some given constants such that Lki,K2{u),u > is a locally 
bounded slowly varying function at 0. Additionally we need to impose a local asymptotic condition on the inverse of 
the transformation z (see below ([H])). 

We state first the result for Pa-x- 
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Proposition 2. Let (X,Y),p G (—1,1) be a bivariate random vector with stochastic representation (1201) . where 
z* : [0, 1] [0, 1] is positive measurable function, R F , and let a G (0, 1] be a given constant. Suppose that 
there exists Op G (l,a/|p|) and for some e G (0,1) the function z[x) := px + z*(x),x G [0,1] is decreasing in 
:— [l/flp — e, l/flp + s] and z{x) < a/ap,\/x G (l/flp, 1]- Suppose that the inverse z^ of z inV^ satisfies 

cd 

(22) z^{alap - d/x) - l/op ^ — 

locally uniformly for d > with c G (0, oo). If further F G GMDAiw) and (|2ip is satisfied with Ki := — l/flp, K2 '■— 
cajop, then Op is unique and 

(23) Pa;. - P{h =h = l}r(7a + I) L k„kA^ / v{x,)) f , 

where a;* := apX,v{x) :— xw{x),x > 0. 

Remarks: (a) If the random variable W appearing in the stochastic representation (PU)) possesses a positive density 
function h continuous at l/op, then under the assumptions of Proposition [2] the asymptotics in (|2ip holds for any 
Ki < K2 with 7a = 1 and 

(24) Lk,,k,{u) - {K2 - Ki)hil/ap), u>0. 

(b) In view of (fT2|l the tail asymptotics oi pa-x given by (|23|) is faster than F(a;). In fact for any constant p. > we 
have 

lim = 0. 

x^oo (w(x))^F(a;) 

Recall that the assumption _F G GMDAiw) implies lima;_j.oo xw{x) = 00. 

c) As it can be seen from the proof of Proposition [2] the local behaviour of z at l/op is crucial. Since we assume that 
z is a decreasing function in 14 the asymptotic of z in (l/op, l/flp + e) is controlled by the asymptotic relation ()22|) . 
Another possibility for z is to assume that it is increasing in (l/flp — e) and decreasing in (l/ctp — e) so that l/flp is 
a locally maximum for z. In this case we can still find the asymptotics of Pa;^;, provided that additionally we assume 
that z~ and z'^ are the inverses of 2; in (l/op — e, l/flp) and (l/op, l/flp + e), respectively such that 

z"(a/ap - l/x) - l/flp , z+ (a/flp - 1/a;) - l/flp ~ — , c_ , c+ G (0, cx)) 

X X 

locally uniformly in a: > 0. 

In order to approximate the joint conditional excess distribution we need an asymptotic approximation as in the 
previous section of Pa.s,ri;x- In the setting of Model B we can approximate another quantity, namely Pa.s,ii,p:x defined 

by 

Pa,s,ri,p-x ■= P{X > x[l + d/v{apx)],Y > ax[l + r]/v{apx)]}, x ^ 00 
with 5,ri £ [0, 00) and Op as above. Note that the difference to Pa.5,ri:x is that above we employ the normalisation 
function v{apx) and not v{x) = xw{x). From the application point of view considering Pa,s.ri.p-x and not Pa,s,ri-x 
is no restriction since the interest is to be able to approximate the joint conditional excess function utilising some 
normalisation function. However, estimating v(apx) leads to complication since also Op need to be estimated. 

Proposition 3. Under the assumptions and notation of Proposition \^if F has a density function f continuous at 
l/ttp, then Op is unique and if further 

cd 

(25) z^{a/ap + d/x) — l/op 

locally uniformly for d > 0, then for any r/, (5 G [0, 00) we have 

(26) Pa,s,7i,p;x ~ -P{^i=^2 = l} ^(ca + l)exp r~r ~? — 7' 

Op V CO + 1 / v(x^,) 

locally uniformly in S,r]. 
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Remark: If we drop the condition ((25|) . then the claim of Proposition |3] stih holds, provided that S > rj > 0. 
We present next two examples. 

Example 4. Let {X, Y) be a bivariate scale mixture random vector with stochastic representation (pO|) where p = 0. 
Consider the function z* given by 

(1 - pe(0,oo),a;e[-l,l]. 

The inverse function of z := z* is z~^{y) = (1 — y^Y^^\ y G [0, 1]. For any a G (0, 1] the equation 

z~^{a/s) = 1/s, s £ (l,oo) 

has the unique solution Up := (1 + 0^)^/^ £ (1, oo). Furthermore ([^ and (^5)) hold with c = a^^^ . 

Let W > with positive density function h being further independent of the positive random variable i? ~ F. If 

F e GMDA{w), then by Proposition [3] 

(27) P.,.„ - P{,,.,..lK-^(l/a,)exp(-i±i^Jl)£|g|. 

Note that if {RIiW, I2RZ* {W)) is a generalised symmetrised Dirichlet random vector, then Ii,l2,R,W are indepen- 
dent and W possesses the density function h{x) — px^~^ g{x^) with g the density function of beta{a,f3). 

Example 5. Under the setup of Example 4, with motivation from the dependence structure of elliptical random 
vectors we redefine z* as 



z*{x) := p^\/\~ x"^, z{x) :— px + z*{x), p,xe{—l,l), p^ := \J \ — p^ . 

First note that z(^p) = 1 is the maximal value of z{x) for any x G [—1, 1]. Hence in order to apply (|26|) necessarily 
a e (p, 1]. It can be easily checked that the assumptions of Proposition [2] are satisfied for Cp := -y/l — 2ap + a? / p*, 
and ((22|) holds with c := (a — p)/(l — ap) £ (0, 00). Note further that ap < a/\p\ and also (|25l) holds. In view of p6| 
we obtain 

28) Pa,5,,,,p;2; ~ P{/l=/2 = l}— 7 T exp 5-^ . 

1 — ap xw(apX) \ Pio.p ' 

In the special case that W"^ ~ 6eia(l/2, 1/2) and P{Ii = 1} = 1/2, i = 1, 2 with /i,/2 independent (the bivariate 
random vector (X, F) is elliptical distributed) we have 

2ap _ 2ap(l - p") 



hil/ap) 



Consequently ([28]) reduces to ([5]) if additionally S = ij = 0. 

4. Three Applications 

Let {X, Y) be a given bivariate random vector. For some high threshold x the approximation of the joint conditional 
excess random vector 

{X^'^Iy^"''^) := {X - x,Y - ax)\X > x,Y > ax, xe{0,oo), a e (0, 1] 

is of some interest in statistical applications if in particular suitable norming constants can be found so that the 
distribution of (^I^l, fI'^"^!) can be approximated by some known distribution function. 

Another interesting problem of the bivariate extreme value theory is the asymptotic independence of X and Y if 
both have distribution functions in some max-domain of attractions. When X and Y are asymptotically independent 
an interesting topic also for application (see e.g., de Haan and Ferreira (2006)) is the estimation of the residual 
dependence index 77. In our last application we give an explicit formula for 77. 

In the light of our findings above we are able to discuss alternative solutions to both problems for the models of 
Section 2 and 3. 
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4.1. Asymptotics of Conditional Excess Distribution. We start by considering the model of Section 2. For any 
s,t positive and some positive scaling function w we have 

P{X["1 > s/w{t), > t/wit)} ^ P^2E^^ x>0. 

Pa-0,0,x 

If X tends to infinity the asymptotics findings of Section 2 to approximate the above ratio. More precisely, under the 
assumptions of Proposition [T] we have 

Pa-,s,t,x Js,t 
Pa;0,0,a; ^0,0 

By the definition Js^t depends on the limit function ^q. Denote by {Ei,E2) a bivariate random vector with positive 
components and survivor function given by Js,t/ Jq.Oj s, i G (0, oo). Then the above asymptotics can be cast into joint 
convergence in distributions, namely if {Xn, n > 1 is a sequence of bivariate random vectors defined in the same 
probability space such that (X„,y„) = (XW, rl"")), n > 1, then we have the convergence in distributions 

(29) (<?(n)xN,g(„)y[-l) 4 iE,,E2), n ^ cx), 

where the scaling function g equals w. 

The limiting random vector has distribution function which clearly depends on a. Further, Ei and E2 can be 
dependent for instance in the setup of Example 3 taking a = 1. In the next model this joint distribution of (Ei, E2) 
is a product distribution which seems to be more relevant for statistical applications. 

Assume next that (X, y), a, Op, p satisfy the assumptions of Proposition [31 For any s^t positive (j26p implies (set 
s := apSjt := apt/a) 

P{X > X + s/w{apx), Y > ax + t/w{apx)} 
P{X > x,Y > ax} 

_ P{X > x{l + s/v{apx}),Y > ax{l + i/v{apx))} 
P{X > x,Y > ax} 



Pay. 

t a„s + cata„la\ , „ n 

exp -^^ TT^ =:exp 

^ ca + 1 ' 



where 



a{ca + 1) 



Consequently, with (X„,y„),n > 1 as defined above (1^^ holds with 

g{x) — w{apx), x > 

and Ei,E2 two independent exponential random variables with mean l/Da^c and l/D* ^, respectively. 
Under the setup of Example 5 

._ 1 - ap , a-p 



ap(l-p2)' ap(l-p2)- 

Thus the convergence in distribution in ([29]) holds in particular if (J7i, ?72) is uniformly distributed on the unit circle of 
K^. We note that the approximation of the conditional excess distribution we do not assume a specific tail asymptotics 
of X and Y. 
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4.2. Asymptotic Independence and MELX-domain of Attraction. A common measure of the asymptotic de- 
pendence between {X, Y) is the tail dependence function 

P{Gi{X)>l-s/x,G2(Y)>l-t/x} 
mm{P{Gi{X) >l-s/x},P{G2iY) >l~t/x}) '- ^^'^*>' 

(when it exists) where Gi, G2 are the distribution functions of X and Y, respectively. If 1(1, 1) — 0, then we say that 
X and Y are asymptotically independent. See for instance de Haan and Ferreira (2006), Reiss and Thomas (2007), 
Hiisler and Li (2009), Das and Resnick (2009), or Peng (2010) for more details concerning modelling of asymptotic 
independence in the context of extreme values. 

We discuss briefly the asymptotic independence for scale mixture distributions with {Ui, U2) specified by our Model A. 
It can be seen by Example 1 that for particular Ui, U2 the limit l{s,t) can be positive, thus asymptotic independence 
does not hold. However, under the setup of Example 2 (|19p implies that l{s,t) = 0,\/s,t e (0,oo), and thus X 
and Y are asymptotically independent and both X and Y have distribution function in the Gumbel max-domain of 
attraction. 

We deal next with Model B assuming that {X,Y) has stochastic representation ([20|) with p G [0,1). The case 
p £ (—1,0) follows with similar arguments. 

In the following we specify the asymptotic behaviour of W and z{W). Explicitly, we assume that for some 71,72 S 
[0,00) 

P{W>1~ s/x} ^ PMW) > 1 - ^ 

™ P{W > 1 - 1/x} ' ™ P{z{W) > 1 - 1/x} ' ^ 

As in Example 1 applying (jl5p we obtain 

P{RW > x} - r(7i + 1)P{W > 1 - l/v{x)}F{x) 

and 

P{Rz{W) > x} - r(72 + l)P{z{W) > 1 - l/v{x)}F{x). 
Next set z{x) :— px + z*{x) < 1, V.t G [—1, 1] and assume that z*{x) < b < l,\/x E [0, 1]. Applying (fT2|) we obtain 
P{Y>x} = P{Ii^ l,l2^ -l}P{R{pW - z*(W)) > x} 

+P{h = -1, 12 = l}P{R{-pW + z*{W)) > x} 
+P{h = l,/2 - l}P{R{pW + z*{W)) > x} 

- P{h = l,h = l}P{Rz{W) > x} 

- P{h = 1, /2 = l}r(72 + l)P{z{W) > 1 - l/v{x)}F{x). 

Similarly, 

P{X>x} - P{/i = l}r(7i + 1)P{VK > 1 - l/w(a;)}F(a;). 

Consequently, in view of (|13p both X and Y have distribution functions in the Gumbel max-domain of attraction 
with the same scaling function w. Let bi{x), i = 1, 2 be defined asymptotically by 

b^{x) ■.= G-\l-l/x), x>l, 

where G~^ is the generalised inverse of Gi, « = 1, 2. In view of (IT2]) we have 



(30) limJ44 . 1. 

x-s-oo 62(2;) 

Furthermore (see e.g., Falk et al. (2004)) 

(31) w{h{x))[G;\l^s/x)~b,{x)] ^ -Ins, Vsg(0,c5o). 

If X,Y are such that the conditions of Proposition [2] hold, then comparing the asymptotics of P{X > bi{x/s),Y > 



b2{x/t)} and P{X > bi{x/s)},P{Y > b2{x/t)} we obtain utilising further ^ and (|30 

l{s,t)^0, Vs,t G (0,(X)). 
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Consequently, X and Y are asymptotically independent with distribution function in the max-domain of attraction 
of a bivariate distribution with unit Gumbel marginals which is a product distribution. 

4.3. Residual Tail Dependence. Modeling of dependence and asymptotic dependence is often done in the frame- 
work of copula, where the marginal distributions are transformed. The asymptotic dependence does not change under 
monotone transformation of marginal distributions. For X, Y with asymptotically independent components it is of 
some interest to quantify the asymptotic independence in terms of some measures. Let Gi,G2 be the distribution 
function of X and Y, respectively. One successful approach to model the asymptotic independence is the estimation of 
the residual dependence index rj G (0, 1] (see Peng (1998,2008,2010), de Haan and Ferreira (2006), Hashorva (2010)). 
So if for some x, y positive 

I^{Gi(X) > 1 — u, G2{Y) > 1 - u\ 

then for any c > 

S{cx,cy) = c'^^'>S{x,y) 

and the function S'„(l,l) is regularly varying with index — 1/77. Other authors refer to rj as the coefficient of tail 
dependence (Ledford and Tawn (1998), Resnick (2007), Reiss and Thomas (2007)). 

As mentioned above in Model A asymptotic independence is not always observed as for instance in the setup of 
Example 1. However, as noted above for Example 2 asymptotic independence is observed. We calculate rj for 
that example. Denote next by G^^^i = 1,2 the generalised inverse of Gi,i = 1,2. Since further lim„_^oo 6i(w) = 
lim„_>.oo 62(1*) = 00, by pip we can write for any y G (0, 1) as w — > 00 

lim '^"(^'^) ^ > 1 - ^/". G2[Y) > 1 - y/u} 

S'„(l,l) P{Gi{X) > 1 - l/u,G2{Y) > 1 - 1/u} 

P{X > G-i(l - x/u), Y > G^\l - y/u)} 



— lim 
= lim 



00 P{X > G^\l - l/u),Y > G^^l - l/u)} 

P{X > biju) -~ lna;/w(6i(u)),y > h2{u) - \^y j w{h2(u))'\ 

P{X > bi{u),Y > b2{u)} 
P{X > bi{u) - \nx/w{bi{u)),Y > bi{u) - lny(l + o(l))/w(6i(u))} 



— P{X >bi{u),Y >bi{u)+o{l)/w{bi{u))} 



U^QO 



Hence by ffT9|) we obtain 

(33) lim §4fTl =■ = vf M + ^)+ + ^^P^"^) 

u^oo bu[l,l) 1(71+ 72 + 1) Jo 

Consequently, since for any c > we have S(cx,cy) — cS{x,y) we conclude that 

r] = 1. 

We consider next Model B. Let therefore {X,Y) be as in our second application satisfying further the assumptions of 
Proposition [3l Since X, Y are asymptotically independent we deal next with the calculation of residual dependence 
index 77. We assume the scaling function w (see the assumptions of our second application) is such that 

(34) lim = c^-\ VcG (0,00). 

u-s-oo 'W(U) 

Since necessarily lim„_>oo uui(m) — 00 we require further that A G [0,cx)). Thus w{x) — x^^^L{x) with L a positive 
slowly varying function at infinity. If A = we assume further that lim„_>oo L(u) = 00. 

Case limsup„^oQ w(u) < cxd: 

In view of ([30| we have for any y G (0, 00) 

b2{u) - lny/w{b2iu)) = 61 (u) + o(l) - (1 + o(l)) In j//w(6i(u)) = 61 (u) - (1 + o(l)) In j//w(6i(u)), u ^ 00. 
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Since further lini„_j.oo bi{u) — oo, (l32t and (p4|) imply for any a;, y G (0, 1) 

lim ■^"(^'^) _ lijn ^{-^ > bi{u) ^ a^-Hnx/w{apbi{u)),Y > bi(u) - a^-Mny(l + o(l))/w(apbi(u))} 
S'„(l,l) ~ P{X > bi{u),Y > bi{u) + o{l)/w{apbi{u))} " 

As in our first application we obtain 

^™ g"n' n = exp{a^-^[Da,c^nx + D* \ny]). 

Consequently 

Since ap > 1, then clearly rj G (0, f]. It is interesting that 77 depends only on ap and A and not on c. 
Case lini„_s.oo w{u) — 00: 

In order to calculate 77 we need to assume further a certain relationship between bi (u) and 62 (u) . In view of (j30p and 
the assumption on w suppose further that 62 is such that 

(35) lim w{b2{u))[b2{u) - bi{u)] ^ ^ G M. 

As above for any a;, y G (0, 00) we obtain as u — > cx) 

Su{x,y) 
lim -— = lim 



S'„(a;,y) _ P{X > bi{u) - a^p-^lnx/wiapbi{u)),Y > bi{u) - {1 + o{l))a^-^{lny + 0/w{apbi{u))} 



«^ooS'„(l,l) «^oo P{X >bi{u),Y >bi{u) + (l + oil))a^-^^/w{apbi{u))} 

= exp{a^p-^[Da,clnx + DlJny]), 

hence again rj — a^^ G (0, 1]. When A = 0, then 77 = 1, otherwise we have rj G (0, 1). Note in passing that neither ^ 
nor c appear in the expression of the residual tail dependence 77. 

For statistical models estimation of 77 is important. In view of our derivation for this model we can estimate 77 by 
estimating first ap and then A. An estimation of A can be obtained as in Hashorva (2010), whereas estimation of ap 
is not as straightforward. In the more specific model of Example 5 ap can be estimated if we estimate p. Estimation 
of Qfp will be discussed in a forthcoming paper. 

5. Proofs 

Proof of Proposition [T] Since R is independent of the bivariate random vector (Ui, U2) and Ui <\ almost surely 
for any (5, 77 G [0, 00) we have 

Pa,5.,rT,x = P{RUi > x{l + S/v{x)), RU2 > ax{l + rj/vix))} 

P{Ui > x/r, U2 > ax/r} dF{r), Vx > 0. 



Let e be a positive constant. The assumption (fTO|) implies that for any constant c G (1, 1 + e) we have P{Ui > 
1/c, U2 > a/c} G (0, 1). In view of 



(36) lim = 0, VaG(l,oo) 

X-+00 F[x) 

and the fact that (Ui, U2) is independent of R we obtain 

(37) pa.5,v-x - / P{Ui> x{l + d/v{x))/r,U2> ax{l + fi/v{x))/r}dF{r) 

Jx{l+S/v(x)) 

(38) = / ^ \is,v{x),S,r,)dF{s/w{x)+x), 



s 



with 

J-/ ^ \ r ^ r,(TT ^ + S/v{x) a{l + r] / V (x)) \ 

f{s,v{x),S,7j) := P\ Ui > ; , t/2 > \ ^ ' ; , s>6. 

I l + s/v(x) l + s/7'(a;) J 
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li S — J] — 0, then pO| implies that 

/(s,l,0,0) = P{C/i > l/(l + s),t/2 >a/(l + s)}, s>0 

is regularly varying at 0. By the max-domain of attraction assumption on F and ([8]) the result for (5 = ry = follows 
easily applying further Potters bounds (see e.g., de Haan and Ferreira (2006)) for the integrand and utilising Lemma 
7.5 and 7.7 in Hashorva (2007). The general case ^ or 77 can be established utilising the result for (5 = = 0, ([9|) and 
the fact that f{s, v{x), S, rf) < f{s, v{x), 0, 0), ^, 77 > 0, and thus the proof is complete. □ 

Proof of Proposition [2] Define next 

x^, := UpX, v{x) — xw{x), Fx{s) :— F{x*[l + s/v{x^,)]), seR, a;>0. 

By the independence of Ii , I2 and RW we may write for any x > 

Pa;x = P{RW > x,pRW - z*{W) > ax}P{Ii = 1J2 = -1} 

+P{RW > x,Rz{W) > ax}P{Ii = l,/2 = 1} 

=: JiP{h = 1, /2 = -1} + J2P{Ii = l,h = 1}. 

If p < 0, then the fact that z* is non-negative implies Pa-x = J2- When p e (0, 1) the assumptions a/p > Up yields 
further 

Ji < P{RW > a/px} < P{R > a/px] = F{a/px). 
Since z(s) < a/up for any s G [l/a^, 1] we have 

P{W > x/r,z{W) > ax/r} < P{W > l/up, z{W) > ax/r} 

< P{W > l/flp, z{W) > a/ftp} 

— 0, Vx, r > 0, x < r < UpX. 

Hence for any £ e (0, 1) 

i<{ap+e)x pca 

J2 = P{W > x/r,z{W) > ax/r}dF{r)+ P{W > x/r, z{W) > ax/r} dF{r). 



J apX J {ap+e)x 

Since (by the assumption) the function z is decreasing and possesses an inverse function z^ in [l/op — e, 1/ap + e] for 
some given e > 0, we have 

f(ap+£)a: 



/ P{W > x/r, z{W) > ax/r} dF{r) 

J anX 



P{W > l/ap[l + s/v{x,)]-\ziW) > a/ap[l + s/v{x,)]-'^} dF^is) 
P{W > l/ap[l + s/vix,)]-\W < z,(a/ap[l + s/vix,)]-^)} dF,(s) 

ev{x^) 

P{cas/{apv{x,)){l + o(l)) >W- 1/ap > -s/iapv{x,)){l + o(l))} dF^(s). 

Hence by the assumptions on W and F applying Potters bound for the integrand and utilising Lemma 7.5 and 7.7 in 
Hashorva (2007) we obtain 

r{ap+s)x 

P{W > x/r, z{W) > ax/r} dF{r) 
°° 'F{x*) 

S'^° exp(-s) dsL_i/ap,ca/ap{^/v{x^)) , , ..^ ■ 

[y\Xi,}} 

In view of p6| Op is necessarily unique, hence applying p2)) as a; — > 00 

Pa;x = F{a/p) + (1 + 0{l))P{h =l2 = l}T{la + l)L ^1/ ap,ca/ apil / v{x,)) j^^^^ + 0(F{{ap + e)x) 

(u(x*))T^" 
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P{Il =h = l}r(7a + l)L^„a,,ca,aS^l^i^*))j^(^^^ 



and thus the resuh foUows. 



□ 



Proof of Proposition [3] By the assumption on the density funetion h W we have that (pij) holds for any 
Ki < K2, Ki, K2 € M with 7^ = 1. As in the proof above for any 5 £ [0, 00) and e > smaU enough we obtain (set 
(car; + J)/(ca+l)) 



Pa.S,, 



{l + o{l))P{h=h = 1} 

P{ca{s - j^)/{apv{x,)){l + o(l)) >W- l/ap > [6 - s) / {apv{x,)){l + o(l))} dF,(s) 



= /2 = 1} /"[(ca + l)s - <5 - 7?] exp(-s) rf /(V«p) n^*) 



X — ?> 00, 



hence the result follows. 



□ 



We conclude this section with the proof of ([18]). 

For all X large and s > by the independence of Si and S2 for any 5i G [0,cxD),i = l,2,s > we may write (set 
G2,x{z) := G2(l - z/x), Si := s ~ S.^, s,x,z € (0, 00)) 

P{Ui > 1 - Si/X, U2>1- S2/x} 

fl 

P{X^Sl > 1 - s,/x -\iv,i = 1, 2} dG2{y) 

o 

P{XA > 1 - s,/x - A;(1 - z/x), i = l,2} dG2Az) 



dG2Az)IG2{l - l/x). 



P{Si > 1 - [s, - X^z]/{xX,),i^ l,2}dG2Az] 
n G.(l - l/x) > ^ " ~ kz]/{xX^,^ = 1, 2} 



Gi(l-l/a;) 

The asymptotic behaviour of Gi , i = 1,2 implies 
P{Ui > 1 - si/x, U2>1- S2/x} 

^ []G,(l-l/x)72 / (max(0,min([5i -Aiz]/Ai,[s2-A2z]/A2)JJ z^^'Uz. 



Now, for (5 = 77 = we may write further 

P{Ui > 1 - s/x, C/2 > 1 - s/x} 



nG.(l-l/a:)72 



i=i 



i(s/Ai,s/A2) 



^min [s — Aiz]/Ai, [s — A22;]/A2 ^ z'^^ ^ dz 



Since A2 > Ai > 0, for any Q < z < s 



s ^ s s — Aiz^s — A2Z 
A2 Ai ' Ai A2 



hence as a; — ?► 00 

P{Ui > 1 - s/x, C/2 > 1 - s/x} 



nG,(l-l/x)72 



X2''' / (s- A2z)''iz'^^-idz + Ap^ / [s ^ Xiz^' z'f^-Uz 



s/Ai 
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